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A. A. Rukhadze 1 '^ and S. P. Sadykova 2 'E3 

1 ' Prokhorov General Physics Institute, Russian Academy of Sciences, Vavilov Str., 38., Moscow, 119991, 
Russia 

2 ' Humboldt- Universitat zu Berlin, Berlin, 124-89, Germany, Newtonstr. 15 
(Dated: 6 January 2012) 

In nonisothermal plasma at temperature T e » Ti diffusion plays decisive role at conditions of smooth 
inhomogeneity when the inhomogeneity size is larger than the Debye radius by more than \jT e jTi times. 
When the inhomogeneity is rather abrupt and the condition is violated, then during the spreading process 
the Maxwcllian relaxation of ion charges becomes significant. Here, we consider these two phenomena together 
and refer to the anomalous character of diffusion, i.e. anomalous diffusion. 

PACS numbers: 52.25.Fi, 52.25.-b 

Keywords: nonisothermal plasma, Maxwellian relaxation, ambipolar diffusion, anomalous diffusion 



I. INTRODUCTION 



II. ANOMALOUS DIFFUSION MODEL 



One can talk of plasma MS IX continous medium when 
its size (and the inhomogenuity size) exceeds signifi- 
cantly the Debye radius. In a case of nonisothermal 
plasma T e » Ti Debye radius is determined through 
the ion temperature and is equal to v<n/wz,i, where 
VTi = \jTijM is the ion thermal velocity and uin = 

Aitqi 2 riio/M is the ion Langmur frequency with qi = 
Ze being the ion charge, riio - equillibrium ion density 
(en e o = qiUio), here e - electron charge. Here and below 
we use the CGS unit system. We write for general case 
but further on we present the results of our calculations 
only for the Hydrogen plasma Z=l. 

The diffusion process is characterized through the 
mean square displacement dependence on time: 



o , s f r 2 n(r, t)df 

< r * >= r A w- ~ * 7 ' 
J n(r, t)dr 



(1) 



If 7 = 1 then diffusion is termed as normal, in all other 
cases one refers to the anomalous character of diffusion: 
7 > 1 - superdiffusion, 7 < 1- subdiffusion. 

When the size of innhomogeneity of ionic component 
perturb ation is larger than the Debye radius by more 
than yjT e /Ti times, then spreading of the ionic compo- 
nent perturbation has a diffusive character and is de- 
fined by the time of the normal ambipolar diffusion 1 . 
If it is less then the Maxwellian relaxation of ionic 
charges becomes dominant during the spreading of 
inhomogenuity^. 

We will introduce a reader below into the theory of 
spreading of perturbation in nonisothermal plasma tak- 
ing into account both mentioned phenomena and show 
the anomalous character of diffusion. 



First of all, we will consider the main equations of the 
studied phenomenon and its constraints. 

The main constraints for the considered below equa- 
tions are: 

bJVi » J 1 ] k V Te » UV e ] V e » k,VTe,W, Vi » kvTi,U) 

. . (2) 

Here v e , Vi arc electron- atom and ion- atom collision 
frequencies, VTe = y/T e /m - electron thermal velocity 
with m being the electron mass. Nonisothermal plasma 
can be only the weakly ionized plasma, in which the 
dominant collisions are charged particle-atom, -molecule 
collisions, u) oc 1/r, where r is the characteristic time 
of spreading and k oc 1/L, where L is the characteristic 
size of spreading. 



1 Equation of motion of electrons and ions has the 
following view: 



M 



0, q e E 



-An e = 0, 



(3) 



Here q e = e, qi = Ze are the electron and ion 
charges, Z < \qt/q\, E are the self-consistent field, 
Vi, v e are ion and electron velocities, Ti = 0, 
T e = const, rii, n e are the ion and electron den- 
sities. 



2 Field equation (Poisson equation) 

divE = A-Ke{n e — Zrii) 

3 Continuity equation of the ion component 

drii 



dt 



+ div riiVi = 



(4) 



(5) 
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4 Equilibrium state satisfying the stationary equa- 
tions Q-jH) at Eq = takes the following view: 



q e n e0 = qifiio = Zen lQ = Const 



(6) 
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Let's consider the small density perturbations : 

n e = n eQ + Sn e , m = n. l0 + 5m, (7) 

where n a Q - denity of the homogenious back- 
ground. These small density perturbations lead 
to the emergence of the weak self-consistent field E. 

System of equations ((3j)-([5|) at © for small perturba- 
tions ([7]) can be deduced to one equation for perturbation 
dm. The rest magnitudes can be expressed through dm- 
Equation for 5m has the following view: 

( CJLi 2 -v s 2 A)^-^ Li 2 ASr H = 0. (8) 

Here uiLi is the equlibrium Langmuir ion frequency, vs 2 = 
ZT e /M is the squared ion sound velocity and = Const. 
Eq. ((8j) represents the partial differential equation of the 
first order with the constant coefficients. 

Let us formulate the Koshi problem for Eq. ([5]), i.e. 
we consider infinite plasma and set the following initial 
condition: 



5m(r,0) = f(f), 



(9) 



and take into account that limr^ooSm = 0. Using the 
Fourier transform 



8m(r,t) = / dkSm(k,t) exp(ifc • r), 



f(k) 



(2tt) 3 



(10) 



dff(r) exp(— ik • r), 



Eq. ([8]) will be reduced to the ordinary differential 
equation of the first time order. 



dSm(k,t) k 2 uj Ll 2 v s 2 jv. 



Ot 



0JLi 2 + k 2 vs 2 



5m(k,t). (11) 



We obtain the solution of Eq. pT|) with the initial con- 
dition 



6m(k,0) = f(k), 



as following: 



uJ Li 2 + k 2 v s 2 J ' 
Making the inverse Fourier transform (|10j) . we can find 

8nAr,t) = [ dkf(k) exp ( ik ■ r — t ^ — % - ' 

J V W ii + k 'VS 

(13) 

Let us analyse the limitting cases of Eq. (|T3"|) . 
A . In a case of rare plasma, when luh — ^ 0: 

w L4 2 



Sm(f,t) = / dkf(k)exp ik - r — t 



= f(r) cxp -t 



(14) 



This solution corresponds to the Maxwellian tem- 
porary relaxation of the charge density perturba- 
tion. 

B . In a case of dense plasma, when lol 2 » k 2 vs 2 ' 

• (15) 



5m(r,t) = / dkf(k) exp I ik ■ r — t 



The integral can be easily solved at Eq. ([9]), when 
f(r) = NioS(r) with the corresponding Fourier 
transform f(k) = p~jg ; where NiQ is the total num- 
ber of perturbed particles at r = 0. The result is 
well knowni: 



Sm(r, t) 



N i0 



8(nDt) 3 / 2 



exp {-r 2 /ADt) , (16) 



where D — vs 2 jv% is the coefficient of ambipolar di- 
fusion and Eq. (| 16[) describes the diffusive spread- 
ing of the localized ion density perturbation (6 - 
functional) . 

Let us determine the < r 2 (t) > Eq. (fTJ) for the 
ambipolar diffusion (|16[) : 

(17) 



< r 2 (t) >= 6Dt, 



here in Eq. ([T]) we used Sm Eq. (TT5|) . As one can 
see this corresponds to the normal diffusion. 

Let us choose the following two initial perturbations /(f), 
the first as representing the infinite wave: 



f(r) = n l0 cos(f ■ fc ) 



(18) 



with the Fourier transform /(fc) = S ^ k k o)+ s ( k + k o) ^ \ 1CTC 

fco is the wave vector; and the second as repesenting the 
localized ion density perurbation mentioned above: 



f(r) = N l0 6(r) 



(19) 



with the Fourier transform f(k) = r^ya ; what will be 
discussed later. We would lik to start with the condition 
®. 

Then, the equation (fT4"]) corresponding to the Maxwell 
relaxation will turn into 

5m(f,t) = mo exp cos(r • k ) 

or in dimcnsionless form 

5m(R,T) ( uj Li T\ - ^ 
= exp cos(i? • e ), (20) 

mo V v % J 

here , fcoro = e*o, R = r/rp, T = tu>Li - dimensionless 
time. 

The equation (|15[) corresponding to the ambipolar dif- 
fusion will turn into 

5m(f,t) = n io exp ( - k ° Vs - ] cos(r • fc ) 
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or in dimcnsionless form 



6ni(R,T) 



= exp 



klDT 

UJLi 



cos(-R-eb), (21) 



here ko is the same wave vector. 

Finally, the equation for anomalous diffusion fp~3|) will 
become 

x r-+\ ( i k 2 Duj Li 2 v s 2 \ - 

dn z (r,t) = n t0 cxp [ -t — 2 , 2 — 2 ) cos ( r ' k °>- 



ULi 2 + k 2 vs' 



or in dimcnsionless form 



8m(R,T) 



nio 



exp 



-t 



k%Duj Li v s 2 



kn 



o v s' 



cos(R-e ). (22) 



As one can easily see in all cases, Eqs. (|2T)|) - (f2"l2"j) . the 
functions get separated by variables T and R. As a result 
we can not define the mean square displacement. This 
means that at the given oscillatory initial perturbaion 
(|18l) no normal diffusion occurs: the initial perturbation 
gets faded down with the time evolution keeping its ini- 
tial functional form in space - cos(/3i?), i.e. only the 
amplitude decreases. This behaviour is peculiar for the 
Maxwellian relaxation ([H]) independent on the type of 
initial condition. However, when we take another ini- 
tial condition (S(r) - functional) for the case of diffusive 
relaxation t|l 5|) then we obtain the normal ambipolar dif- 
fusion ([To]) . 

Finally, let us find the solution for the anomalous dif- 
fusion Eq. (|13[) with the second initial condition Eq. 
(|19|) as representing the localized ion charge perturba- 
tion. Having solved this equation and having made some 
simplifications, we have got the following: 



6nj(r,t) _ r 3 D 
N i0 /r 3 D ~ (2tt)3 



dkexp [ik-r — t 



1 D 



dkf{k)k 



k 2 u Ll 2 vs 2 /vi 
UJLi 2 + k 2 vs 2 
sin(fcr) 



x exp —t 



k 2 v 2 u> 2 Li /vi 



k 2 v 2 



+ H^) r r> CX P ( — * 



^expl-t- 



dkk 



sin(fcr) 



(23) 



III. RESULTS AND DISCUSSIONS 

In the Figure [1] a) the 3 dimensional view of the rela- 
tive value of density perturbation in dependence on time 
T and absolute value of the radius vector R is shown. 
Namely, Fig. [T] a) represents the anomalous relaxation 
as the combined effect Eq. (|2"2"]l of the Maxwellian re- 
laxation Eq. (|20p and ambipolar diffusive relaxation Eq. 
(|2ip . In Fig. Q] b) the comparison among the mentioned 
relaxations and the localized ion density perturbation re- 
laxation Eq. (|16|) at t = 1/uiLi is shown. 



As one can easily see the behaviour of the Maxwellian 
relaxation Eq. (|20[) . ambipolar diffusive relaxation Eq. 
(|21j) . and the anomalous diffusion Eq. (fT6|) is oscillatory 
and fades down with the time as the exp(— at) cos(/3i?) 
function. In a case for the localized ion density pertur- 
bation relaxation Eq. (|16[) the perturbations fades down 
as the cxp(-jR 2 /t)/t 3 /' 2 . 

At the given oscillatory initial perturbation (TT5)) no 
normal diffusion occurs: the initial perturbation gets 
faded down with the time evolution keeping its initial 
spatial functional form - cos(/?i£), i.e. only the ampli- 
tude decreases, and in this way revealing the anomalous 
nature of the diffusion in nonisothermal plasma. This 
behaviour is peculiar for the Maxwellian relaxation (TH|) 
independent on the type of initial condition. However, 
when we take another initial condition (5{r) - functional) 
for the case of diffusive relaxation (|15p then we obtain the 
normal ambipolar diffusion (|16l) . 

In the Figures [31 L3 HI an d E] comparison between 
the obtained relative density perturbation relaxation for 
anomalous diffusion Eq. (|23[) with S(r) as initial per- 
turbation (|19[) and normal ambipolar diffusion Eq. (|16[) 
and corresponding variances at the different ion Lang- 
muir frequencies and different time moments are shown. 
In the Figs. [H [3] the difference between the obtained 
curves for anomalous diffusion and the normal ambiplar 
diffusion is quite significant, whereas with an increase of 
either ion Langmuir frequency or time the anomalous dif- 
fusion gets to turn into the normal ambipolar diffusion. 
This is also demonstrated in the corrcspoding Figs, d) 
where dependence of the variance on time and the fitting 
curves for the corresponding different Langmuir frequen- 
cies is shown. The variance was estimated using Eq. (TTJ) . 
As we can see at the lowest frequency there is almost 
no diffusion Fig. [2] d), the variance does not depend 
on time. However, with further increase of the ion fre- 
quency the diffusion becomes anomalous < r 2 (t) >^ at 
(a=Const), whereas when the frequency and time mo- 
ments increase the diffusion starts to turn into the nor- 
mal: < r 2 (t) >= at (fT7|) . 



ACKNOWLEDGMENTS 

One of the authors S.P. Sadykova would like to thank 
her coauthor A. A. Rukhadze for his fruitful contribution 
to the work and discussions on the presented topic. S.P. 
Sadykova would like to express her gratitude to her father 
P. S. Sadykov for his financial support of the work and 
for being all the way the great moral support for her. 

X A.F. Alcxandrov, L.S. Bogdankcvich, A. A. Rukhadze, in Prici- 
ples of Plasma Elektrodynamics, Springer Series in Electro- 
physics, Volume 9, edited by G. Ecker (Ruhr-Universitt Bochum) 
(Springer- Vcrlag, Berlin, Heidelberg, New York, Tokio, 1984). 

2 M.V. Kuzelev A. A, Rukhadze, in Basic of Plasma Free Elektron 
Lasers (Edition frontier, Paris, 1995). 

3 M.V. Kuzelev A. A, Rukhadze, in Methods of wave theory in dis- 
persive media (World Scientic Publishing, Singapore, 2010) 



4 




FIG. 1. a) The anomalous diffusion Eq. 1)221) and b) comparison among the anomalous, Maxwellian relaxation Eq. (|20l) . 
ambipolar diffusive relaxation Eq. (|21[) and the localized ion density perturbation relaxation Eq. (|16[) at t = 1/ujlx is shown. 
Here, eo\\R, T e = 10 eV, n i0 = 10 12 cm" 3 , u L i = 1.3-10 9 s _1 , Vi = 1.2-10 9 s _1 , n = 3-10 16 cm" 3 (atom density), r D = 2.3-10" 5 m. 



4 L.D. Landau, E.M. Lufshitz. in Hydrodynamics (Nauka, Moscow, 
1982) 



5 





FIG. 2. Comparison between the obtained the relative density perturbation relaxation for anomalous diffusion Eq. (|23p and 
normal ambipolar diffusion Eq. (|16[) at the different time moments T = tVi but fixed Langmuir frequency uiLi = 1.6 ■ 10 7 s _1 , 
riia — 1.5 ■ 10 s cm -3 and its variance Eq. fl]) together with the variance fit d) are shown. Here T e = 10 eV, Vi = 1.2 ■ 10 9 s _1 , 
no = 3 • 10 16 cm -3 (atom density). 
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- Fitting, <FT>=0,641 16-0,24226 T+0,25896 f 

- Anomalous diff. 



3 

T=tv 



C) 



d) 



FIG. 3. Comparison between the obtained the relative density perturbation relaxation for anomalous diffusion Eq. (|23p and 
normal ambipolar diffusion Eq. (|f 6[l at the different time moments T = tfi but fixed Langmuir frequency uiLi = 1.3 • I0 9 s~ 1 , 
mo = f0 12 cm -3 and its variance Eq. JT]) together with the variance fit d) are shown. The Hydrogen plasma parameters are 
the same as in the Fig. [2] 
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FIG. 4. Comparison between the obtained the relative density perturbation relaxation for anomalous diffusion Eq. (|23p and 
normal ambipolar diffusion Eq. (|f 6[) at the different time moments T = tvi but fixed Langmuir frequency uiLi = 3.3 ■ fO 9 s _1 , 
mo = 6.25 • fO 12 cm -3 and its variance Eq. |T} together with the variance fit d) are shown. The Hydrogen plasma parameters 
are the same as in the Fig. [2] 
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FIG. 5. Comparison between the obtained the relative density perturbation relaxation for anomalous diffusion Eq. (|23|) and 
normal ambipolar diffusion Eq. (|f 6[) at the different time moments T = tvi but fixed Langmuir frequency uiLi = 9.7 ■ fO 9 s _1 , 
mo = 5.6 • fO 13 cm -3 and its variance Eq. ([TJ together with the variance fit d) are shown. The Hydrogen plasma parameters 
are the same as in the Fig. [2] 



